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In Physica D 91, 223 (1996) [H, results were obtained regarding the tangent bilurcation of 
the band edge modes {q = 0, n) of nonlinear Hamiltonian lattices made of A'^ coupled oscillators. 
Introducing the concept of partial isochronism which characterises the way the frequency of a mode, 
ui, depends on its energy, e, we generalize these results and show how the bifurcation energies of 
these modes are intimately connected to their degree of isochronism. In particular we prove that in 
a lattice of coupled purely isochronous oscillators {w{e) strictly constant), the in-phase mode {q = 0) 
never undergoes a tangent bifurcation whereas the out-of-phase mode {q = tt) does, provided the 
strength of the nonlinearity in the coupling is sufficient. We derive a discrete nonlinear Schrodinger 
' equation governing the slow modulations of small-amplitude band edge modes and show that its 

' nonlinear exponent is proportional to the degree of isochronism of the corresponding orbits. This 

^ , equation may be seen as a link between the tangent bifurcation of band edge modes and the possible 

O ■ emergence of localized modes such as discrete breathers. 
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I. INTRODUCTION 



Since almost fifteen years now, properties of discrete breathers (DB) in nonlinear translationally invariant Hamil- 
L^j, tonian lattices are under intense investigation. General features regarding these time periodic, spatially localised 
excitations are well understood and have been documented in several reviews 0, 0, 0; S IE IE H- Mathematical 
, proofs of their existence go back to a paper by R.S. MacKay and S. Aubry |E| which considers lattices of interacting 
^ oscillators. It shows the possibility to continue single site oscillations in the decoupled (so-called anti-continuous) 

limit to nonzero coupling between the oscillators provided the corresponding orbit stays out of resonance with the low 
• amplitude lattice modes. This, in particular, requires the oscillator (onsite) potential to be non isochronous, that is, 
' to possess orbits whose frequency varies with the energy. Other works, either generalising this method 9], or using 
T— I different approaches JO, 11, 12] have added to the variety of rigorous existence proofs of discrete breathers. 

One of us has performed an analytical study of the way band edge modes (BEMs) may bifurcate (tangently) to give 
, rise to new periodic orbits breaking the translation invariance of the lattice |13| . It is a common conjecture indeed 
(~| ' that discrete breathers are among the orbits bifurcating from these plane waves. The bifurcation analysis investigates 
' ^ ' the possible existence of almost extended discrete breathers and, in this respect, is complementary to MacKay's and 
^ Aubry's theorem which proves the existence of strongly peaked ones in networks of weakly coupled non isochronous 
oscillators. It provides the critical energy at which tangent bifurcations of BEMs possibly occur according to the 
•rH , mode under consideration (see also But it is restricted to the generic case where plane wave orbits do not bear 

' any degree of isochronism. 

In order to discuss further the specific properties of partially isochronous BEMs, we now define this concept more 
precisely. Let us first remark that in one-dimensional (ID) convex potentials, the motion is always periodic. To any 
given energy corresponds a unique orbit whose frequency is determined by the features of the potential (basically, its 
shape) . We will say an orbit to be isochronous if its frequency does not depend on its energy. Given the one-to-one 
correspondence between the potential and its orbits, if the motion is isochronous the potential can be said to be 
isochronous as well. 

In higher dimension however, the concept of isochronous potentials becomes ambiguous due to the possible presence 
of several families of periodic orbits, each of which having its own energy-frequency dependence. Isochronism is then 
a property of a particular family of orbits rather than a property of the potential itself. Let us notice though, that 
to give birth to a family of isochronous orbits, the potential must fulfill certain conditions. These constraints simply 
vary according to the motion under consideration, and may be different for different families of periodic orbits of one 
and the same potential. 

The most famous example of a ID isochronous potential is the harmonic well V{x) = whose frequency uj is 

well known to be energy- independent. Nevertheless, isochronism is not the privilege of the latter and it can be shown 
that appropriate shears of the parabolic curve produce other non .symmetric isochronous potentials V{x) ^ Vi^—x) 
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Now, for generic convex potentials, the frequency of a given periodic orbit can be expanded at low energies E (bottom 
of the potential) as a power series in E. Its behaviour is generally linear with E around the equilibrium position. We 
will call "partially isochronous" or, more precisely isochronous up to order n, orbits whose frequency behaves instead 
as a nonlinear function of the energy when expanded around E = 0. Typically, U!'^{E) = ujq + 'fnE^ + o[E)^, n > 2, 
ujQ > and 7„ ^ 0. For n = 1, we recover the case of non isochronous motions or, equivalently in our terminology, of 
orbits isochronous up to order 1. Completely isochronous orbits verify lj'^{E) = Wq. 

The aim of the present paper is to extend the results 13] to the class of potentials rendering certain periodic 
motions partially or completely isochronous. As we will see, it turns out that the tangent bifurcation of band edge 
plane waves is intimately related to the low-energy behaviour of their frequency (hence the above definitions). This 
remark will enable us to generalise the perturbative analysis performed in [T^ and to express the bifurcation energy 
of BEMs in a very simple way valid for any potential. An immediate conclusion of this generalisation will be that, 
families of discrete breathers bifurcating from partially isochronous BEMs possess energy thresholds, even in ID. This, 
to some extent, shall be confirmed by a "multiple-scale" analysis which leads to a discrete nonlinear Schrodinger 
equation (DNLS) whose nonlinear exponent is related to the degree of isochronism of the BEM. The existence of 
energy thresholds of DBs for such DNLS models has been confirmed J^] and has been proved by M.I. Weinstein [l^ . 

The paper is organised as follows: in section m we present a low-energy perturbative solution of the in-phase 
and out-of-phase modes and in section IIIII wc use them to derive the constraints imposed to ID or 2D potentials 
to enforce a partial isochronism of these motions. In section IIVI we study the linear stability of these orbits and 
derive an expression for the energy at which they undergo a tangent bifurcation. Finally, in section^ we discuss the 
implications of these results and present a brief multiple-scale analysis which corroborates our conclusions. 



II. EQUATIONS OF MOTION 
A. Generalities 



All along this paper, we will investigate the dynamical properties of a lattice described by the following Hamiltonian 



N 



+ V{Xn) + W{Xn+l - Xn) 



(1) 



with periodic boundary conditions Xn+N — a^n- For the sake of simplicity, we consider an even number of sites N 
ranging from 2 to infinity. The onsite (V{x)) and the interaction {W{x)) potentials are both assumed to possess a 
minimum at a: = around which they can be expanded as 



OO -J OO 

V{x)^y2-v,x'^ ; W{x)^y2-<j>^xf^ 



(2) 



fj.=2 



11=2 



The first coefficients of these expansions, V2 and (/)2, represent harmonic frequencies and are assumed to be strictly 
positive. The Hamiltonian equations of motion for are given by 



•^n Pn 1 

Pn = -V'{Xn) -W'{Xn- Xn-l) +W'{Xn+l- Xn)- 



Let us introduce the normal coordinates 



1 ^on 27rZ , , Af ^ A^, 



Their properties are 



iq+2-n: 



and inverting the transform yields 



}q and g_<, = 0* {xn e M), 



(3) 

(4) 

(5) 
(6) 
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Rewritten in terms of normal coordinates, equations Q read now 
where 



1 ^ 
N ^ 





A linearization of Fq{Q) around Qq — leads to the equations of motion of a harmonic lattice, namely 



where 



V2 + 402 sin' 



(I) 



(7) 



(8) 



(9) 



(10) 



represents the squared frequency of each linear mode q (hence the additional subscript 0, the frequency of the nonlinear 
mode being denoted by LUq). 

In what follows, we will be interested in the stability of two particular nonlinear modes corresponding to the natural 
continuation of the linear q — and q = n modes defined by © • These nonlinear modes are periodic solutions of (0) 
which converge to their respective linear modes as their energy tends to zero. Notice that the linear frequency luq.o 
of the in-phase mode is always nondegenerate and because the number of sites N is even, the linear frequency 10,^,0 
of the out-of-phase mode is nondegenerate as well. All other modes q 0,tt are twofold degenerate (wg,o = ^-q,o)- 
In the next two sections, we define more precisely the two nonlinear in-phase and out-of-phase modes and evaluate 
them perturbatively by means of a Poincare-Lindstedt expansion carried out at low energy. 



B. In-phase mode (orbit I) 

1. Equation of motion 



Oscillators are said to be in phase when they perform identical periodic motions. This corresponds to 

Qg^QoSq,o (11) 

where Sq^q' = 1 if g = [2tt] and else. The previous expression is a solution of the equations of motion 10 provided 

Qo + V'iQo) = 0. (12) 

The solution Qq = QoSq o represents the in-phase periodic orbit. We call it orbit I. The total energy of the lattice 
evolving according to orbit I is 

Ej - H{{xn - Qo}) - N ^ v{Qo)^ . (13) 

We will use an energy density (or energy per site) rather than the total energy Ej to describe this orbit. It is given 
by Sj = Ej/N and represents the energy of the oscillator Qq evolving according to H12|l . 



2. Solution at low energy 

Eq. H12|) represents the motion of a single oscillator in the potential V{x). According to our assumptions regarding 
the latter, the potential is convex in a; = and at small energy, the motion is bounded and thus periodic. We can 
solve for it in perturbation by expanding the solution as a Poincare-Lindstedt series (see e.g. H^)- For this purpose, 
let us first define a new dimensionless time r = ui{e)t, where Lj{e) is the frequency of Qq as a function of its energy 
e (for the sake of clarity we have dropped the subscript I). The corresponding period is T(e) = 27: /Lj{e). We show in 
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Appendix 1X1 how this period can be expanded as a power series in energy. Let us also define X = (2e /v2Y^'^ ^ as well 
as the new dimensionless quantities T{X) = T (e) y/v2 / {"^tt) and V{x) = V{x)/v2. Equation II12|) now reads 

^^^^ + f\X)V' (Qo(X, r)) = 0, (14) 

where we have explicitly mentioned the energy dependence of Qq on X. To solve this last equation, we expand 
Qq{X, t) as a series in X around 0: 

oo 

Qo(^,r) = 5]Q("HT)X". (15) 

n=l 

We choose without lost of generality an initial condition such that Qo(-'^jO) = which means that Q[){X,{)) is a 
turning point of the potential V defined by the relation V{Qi){X,Q)) = e V{Qi){X,Q)) = X"^ /2. Inverting the 
previous relation gives 

oo 

Qo(^,0) = X + ^(T„X", (16) 

where the first odd coefficients (t„ are given in appendix^ According to the relation above, the initial conditions for 
the functions Qo"''(t) are 

Q^'^(O) = 1 ; Q^"^(0) = a„ , Vn > 2 and q(")(0) = , Vn > 1. (17) 
Moreover, according to l|A6|l . we know the explicit form of T{X) which reads, 

f{X)^\ + Y.^2kX^'' where f2k=<T2k+i^^^. (18) 
fe=i ^ '" 

Reinstating H18|) and H15|l in (|14|l and expanding as a series of X, we derive a set of differential equations for Q^q\t). 
This way is clearly related to the Poincare-Lindstedt method except that the preliminary calculation H18|) of the period 
T{X) automatically removes all secular terms from equation H14|) . 

Using the results of Appendix^ we may derive for example the first differential equations involving Qq^"* (t), Qq^'' (t) 
and Qfir). 

Q'^^\t)+q'-^\t)^Q, (19) 
(r) + g(^) (r) + a2 [g^'^ (r)] ' = , (20) 

Q'i^ (r) + Q^o^ {r)+{^-al- Q^'^ (r) + as [q^^^ (r)] ' + 2a2Q^'^ (^)q(2) (^) ^ q , (21) 

where the double dot stands now for a differentiation with respect to t and a„ = Vn+i/v2- Solving this system 
together with the initial conditions Ijl?!) . we obtain 

Q^^\t) = cos(t), (22) 
Q^^\t) = (cos(2t) ~ 3) , (23) 

Qo\r) = ^ {2al + 3as) cos(3r) + (^-^al - ^ag) cos(r) . (24) 

where we have used a2 — —a^j'i and the value for uj, given in Appendix 1X1 

The first term of the expansion (|15|l . Q^q\t) ~ cos(r), represents the linear (harmonic) part of the in-phase motion. 
Higher order corrections in X", n > 1 (i.e. et"/^)) stem from the nonlinearity of the onsite potential. It is important 
to notice that they appear also for isochronous onsite potentials. In other words, isochronous motions are stricto 
sensu anharmonic. 

Reinstating the results for the Qq"'' (t) in (|15|l provides the general perturbative expression for the continuation of 
the linear in-phase motion, that is, a perturbative expression for the nonlinear in-phase mode. 
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C. Out-of-phase mode (orbit II) 

1. Equations of motion 

If the potential V{x) is not symmetric {V{—x) 7^ V{x)), its Taylor expansion around contains at least one nonzero 
odd cocfBcient. This has no influence on the previous result concerning the in-phase motion because any oscillator 
of the chain performs the same motion in the same time. This reduces the set of N equations Q to a single one 
(|12|l . representing the equation of motion of a single oscillator in the onsite potential V. But as soon as we are 
interested in an out-of-phase like motion, we have to consider a dimerization of the chain, each dimer being made 
of two neighbouring units oscillating in opposite phase. The lack of symmetry of V induces two different motions 
to the right and to the left. This prevents us from finding a pure out-of-phase solution to which would imply 
Qq — Q-nSq^-K Or lu real space X2n — —X2n+i- Instead, we can look for a solution of the type 

Qq = QoSqfi + QTrSq^T, (25) 

involving both in- and out-of-phase variables, the others being zero. Using ©, we obtain Xn — Qo + (— l)"Q7r or 
X2n ~ Qo + Q-K and X2n+i = Qa — Q-n- Adding and substracting the equations of motion for X2n and X2n+i, wc finally 
get 

Qo + ^ [V'iQo + Q.) + V'iQo - Q^)] = , 

Q^ + l [V'iQo + Q.) - V'iQo - Q^)] + W'{2Q^) - W'{^2Q^) = 0. (26) 
The total energy of the system evolving according to orbit II is 

Ell = H{{x2n = Qo + QiT, X2n+l = Qo — Qtt}) 

= ^[qI + QI + y{Qo + Q^) + ^(Qo - Q.) + W{2Q^) + W{-2Q^)) . (27) 



2. Solution at low energy 



System represents a dimer. Applying the method of the first section, we are able to derive its time-periodic 
solution at low energy by requiring that the corresponding orbit converges towards the linear out-of-phase mode as 
the energy vanishes. However, at variance with the in-phase mode, we cannot provide an explicit expression for its 
period in terms of the Taylor coefficients of V and W . Removing the secular terms from (|26f) yields simultaneously 
the expressions for the period and for the motions Qo{t) and Qtt^t). 

Contrary to the previous case, we won't use an energy expansion for the diverse quantities to be calculated but 
merely a small amplitude expansion. The small amplitude is denoted by y. We define a dimensionless time r = uj{y)t 
where uj{y) stands for the frequency of orbit II as a function of its amplitude y. The corresponding period is denoted 
by T{y). Moreover, we define the two dimensionless potentials 



V{x) 



V{x) 



V2 + 402 



oc 

EOr. 
i 

n=2 



W{x) 



W{x) 

V2 + 402 



1 

n=2 



(28) 



as well as a dimensionless period T{y) = ^Jv2 + 402 /2t:. Notice the slightly different definition of the coefficients 
oin = Vn+i/{v2 + 402) of this scctiou as compared to the previous one. We nevertheless keep the same notation for 
the sake of clarity given that these coefficients are still related to the onsite potential V{x). Functions Qq, and 
T{y) are expanded as follows. 



Qo{y,r)^J2Qo'\r)y'' ; Q.(y,r) = ^ Qt")(r)2/" ; f (y) = ^ f; ; 



(29) 



Equations of motion H2t)|l now read 



f'{y) 



V'iQo + Qrr)+V'{Qo~Q.) 



Qr.+T^{y)\i: V'iQo + Q.)-V'iQo-Q.) +W'i2Q^)-W'i-2Q^) \ =0 



(30) 
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The double dot denotes the derivative with respect to r. Inserting H29() in the previous system and collecting the 
terms of same order in y gives rise to a set of differential equations involving the functions Qg"' (r) , (r) as well 
as the unknowns Tn (to be determined by requiring the removal of secular terms). Solving them, it is not difficult to 
obtain the following general features for the motion: By choosing a proper origin of time the solution can be made 
time-reversal symmetric (the only nonzero Fourier coefficients of Qq"^(t) and (5i"''(r) are even). Moreover, (5o"^(t) 
and qI"^ (t) are respectively even and odd in y. Finally, Tijj) is also even in y. Thus, we have 



OO 

n=l 



OO 



^'\r)y 



2n+l 



T{y) 



E^^ 



2ny 



(31) 



where 



gW(r)=cos(r), 



1 



)(3) 



To 



2 
1 

32 

1 (3ai - 



Ol2 



cos(2t) 



4- 

("3 - 



8)«1 



1 

ai 



■P3 



2al 



4 (ai - A)ai 8 



ai ~ 4 
(a3 + /93). 



cos(3r). 



(32) 



to give a few. Notice that the even coefficients (32n are absent from the equations of motion (|30|l and consequently 
from the expressions (|32|l . (3i has been eliminated thanks to the relation /3i +ai — 1, which follows from the definition 
of these coefficients in terms of V2 and 02 (see eqs. and (EHJ)- 

Reinstating properly the expressions for Qo{y, r) and Qj^iy, t) in (|77|l allows us to derive a perturbative expansion 
for the energy density e = En/N of the nonlinear out-of-phase mode in terms of its amplitude y. As for the in-phase 
mode, it is convenient to define a quantity Y — {2e/{v2 + 4:(f>2)y^'^ playing the same role as X in the previous section. 
We then obtain 



2e 



V2 + ■ 



9 , 



-P3 



68ai + 96)a: 



ai{ai - 4)2 



y 



0{y'). 



(33) 



The first term of this expansion corresponds to the harmonic limit. 



III. ISOCHRONISM 



In the previous section, we have described the way to obtain a perturbative energy expansion for both the in- and 
out-of-phase motions, as well as for their respective periods. We are thus in a position to express the conditions 
required for a (partial) isochronism of these modes. For an isochronism of order n, this is easily achieved by cancelling 
all coefficients of the energy expansion for the period up to order n — 1. 



A. In-phase mode 



Due to its integral representation, the period of the in-phase mode can be explicitly expanded as a power series 
in energy (see Appendix ^) . Obtaining an in-phase mode isochronous up to order n then amounts to zeroing the 
coefficients <72k+i, 1 < fc < n — 1. This has been done in Appendix up to order n = 4. The set of equations thus 
derived induces some constraints on the Taylor (vk or ak-i) coefficients of the potential V. It leaves nevertheless 
an entire freedom on the choice of odd coefficients (w2fc+i, fc > 1) as already remarked in For a (ID) potential 
isochronous up to order n, the even coefficients V2k are then determined by the odd ones V2k+i up to fc = 2n. The 
rest of the expansion is free. 

To illustrate this, let us derive the most general expansion of a {C°°) ID potential isochronous up to order 4. Using 
the results of appendix 1X1 and V{x) — uj'^V{x), we obtain: 

1 2 "2 3 5 2 4 "4 5 1 / 56 4 14 \ g 



Vix) = -x' + —x^ + —aix^ + —x'' + -l-—al + —aia2]x 



"6 7 1 /24 592 3 36 , 848 o ^ ak-i . 

' ' -"6^2 - -^a^al + —a\ + —a\ J + E ^ (34) 



A;=9 
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This expression is equivalent [4J| to Eq. (34) of |23 |. 

B. Out-of-phase mode 

As for the out-of-phase motion, no cxpUcit (integral) expression is available for the period. However, we have shown 
in the previous section how to calculate it as a function of the small amplitude y by removing the secular terms at 
each step of the calculation. The period could of course be converted as a series in energy rather than expressed as 
a series in amplitude. Such a transformation would be pointless however. Indeed, even if an isochronism of order 
n has been previously defined through T = Tq -I- r2„e" -f o(£") as the energy tends to zero, the dependence of the 
amplitude on energy H33|l shows that this statement is equivalent to T = To -|- T2n[{v2 + 402)/2]"y^" -I- o{y^"'). It is 
then equivalent to zero the coefRcients of the y or e expansions up to order n — 1 as it produces the same constraints 
on the Taylor coefficients of the potentials V and W, although these two series have different coefficients. 

Let us finally provide the reader with an example to illustrate the method described in the previous section. If we 
require the out-of-phase mode to be isochronous up to order 3, we obtain the following contraints on the coefficients 
a and f3 (see eqs. (|28|l for their definition): 



2 



/?3 = -"3 



2 (3ai - 8)a 

3 (ai — 4)q!i 

2 (96 + 15a?-56ai)^4 9 (32 + 5a? - 24ai) ^ ^ 6 (5ai - 12) ^ 



f35 = -a5 + — 37 TT^ "2 --7 3? Tvz "i"2"3 + --37 — aia2a4 (35) 

15 a^(Q!i— 4)^ 5 Q;j(ai— 4)^ 5Q;^(ai— 4) 

Notice that for the sake of simplicity, we have kept the same notation at for the coefRcients of V for both modes. 
However, the definition of these coefficients is mode dependent as they represent Vk+i/v2 for the in-phase mode and 
Z)fe+i/(w2 -l-4<^2) for the out-of-phase mode. As already displayed in the above relations, it can be shown in general that 
the isochronism of orbit II is easily expressed through a relation of the type /32n+i = /({"fe}) {k £ {1, . . . , 2n -I- 1}). 
This means that once the coefficients {ak} have been chosen for the onsite potential V, the isochronism of the out- 
of-phase mode fixes the even part of the interaction potential. The odd part of W remains completely free as it does 
not enter the equations of motion. 

IV. TANGENT BIFURCATION OF ORBIT I AND II 

A. Statement of the problem 

The system of equations ^ is of the form Q + F{Q) — 0, where Q and F{Q) denote two vectors of components 
Qq and Fq{Q) respectively, q £ {0, 2tt/N, . . . , 2(A^ — l)7r/A^}. A perturbation yy of the system around the solution Q 
gives rise to the following variational system 



To evaluate the Jacobian matrix, we use ijHJl and find 



?7 + DF{Q)t] = (36) 

dFg 
dQk 



where DF{Q) is the Jacobian matrix of F evaluated in Q whose components are DFqk{Q) — §^iQ)- 



n— 1 



.(e-'fc - 1)W" I ^(e-'« - l)e-"?'"Q, 



(37) 



We will now use the method already described in in order to obtain the energy thresholds, if any, at which 
orbits I and II become unstable and bifurcate tangently to give rise to other types of periodic orbits which break the 
translational invariance of the lattice. Before to proceed, let us briefly recall how the method works. 

Once the periodic solution for Q has been introduced in the Jacobian matrix DF{Q), the variational system H36|l 
presents itself as a vectorial Hill's equation for the perturbation rj. This type of systems is known as parametrically 
excited as the Jacobian matrix generally depends on several parameters. In our case, such parameters are the energy 



8 



(or the amplitude) of the solution Q as well as the frequencies of the modes we are interested in. Once expanded as 
a Fourier series, the Jacobian matrix may be decomposed into a static (dc-) part (its zero mode) and a driving (ac-) 
part. 

A paradigmatic example of Hill's equation is the Mathieu equation x{t) + (5 + 2e cos(2t))a;(i) — 0. For a compre- 
hensive treatment of this equation the reader is invited to consult ;21j, for example. In one dimension (the Jacobian 
matrix is reduced to a single element in this case), the parameter 5 plays the role of the static part and 2ecos(2t) the 
role of the driving. It is known from the stability analysis of this equation that the behaviour of its solution varies 
according to the values of the parameters 5 and e. The solution can be either stable, unstable or periodic. In the 
ei5-plane [5 as a:- axis and e as y-axis), the regions of instability present themselves as tongues (the so-called Arnold's 
tongues) starting from the J-axis at the values (5„ = n^, n g N and widening as e increases. In these regions, the 
motion is unbounded whereas outside, it is stable (bounded). 

Of particular interest are the boundaries of such regions called transition curves that separate stable from unstable 
motions. Along this curves, the solution is periodic of period tt (n even) or 27r (n odd). For small e values, a 
perturbative treatment of the Mathieu equation which consists in expanding both x{t) and the parameter 5 as series 

in e allows for the determination of the transition curves of the form 5n — rP' + S;^ni°^'- coefhcients A^^^" 

depend on the tongue (n) as well as on the branch (that is, the boundary) we are interested in. It can be shown 
that one of these branches is related to time reversal symmetric solutions r] (denoted by the subscript s) whereas the 
second one is associated with time reversal antisymmetric solutions (a). 

Let us suppose now that wc fix the value of 6 close to a transition point 5n — ■ At e = 0, the point corresponding 
to the state of the system in the parameter space is located in a region of stability. Let us increase the value of e at 
fixed 5. If the corresponding vertical line crosses the transition curve nearby, the solution x{t) becomes unstable above 
the crossing point. And right at the intersecting point, the solution is periodic. Although the variational equation 
is more complex than the Mathieu equation, we will proceed along the same lines as those described above to 
derive the energies at which the modes (in- or-out-of-phase) become unstable and bifurcate. 



B. In-phase mode 



1. Stability analysis. 



Evaluated along orbit I, the Jacobian matrix is diagonal 



y"(Qo)+402 sin^ 



All perturbations decouple from each other and their equations of motion are 

^9 



^^"(Qo) + 4(/.2sin2(|)l,;, = 0. 



(38) 



(39) 



As already stated in , the perturbation 770 describes the continuation of orbit I along itself. It cannot be responsible 
for a bifurcation of Qq as it simply operates a shift in time or modifies the energy (or the frequency) along the one- 
parameter family. We then look for the perturbation able to give rise to the required tangent bifurcation. The first 
to occur will be for the closest (linear) frequency to the linear in-phase frequency, that is for Qc = 27t/N. 

Let us now rewrite the system made of the dimensionless equations of motion for Qo{t) and rjq^ir). We obtain 



+ f\X) \V" (Qo(^,t)) - aI 77g^(X,r) = 0, 



where 



"2 . 2 
— sm 



V2 



^0.0 



V2 



(40) 
(41) 

(42) 



The notations used above are the same as in section lTlBI The energy dependence of all quantities has been emphasized 
through X. 
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2. Arnold tongue corresponding to a tangent bifurcation. 

Recall that due to the thiie scaling r = uj{e)t the period of Qo{X,t) is now 2tt. So that looking for a tangent 
bifurcation of this orbit through riq^{X,T) requires r]g^{X,T + 27r) = riq^{X,T) as well. However, this condition in 
itself, although necessary, is not sufficient for our purpose. We have indeed to require the frequency of rjq^ {X, r) to 
be the same as (and not a multiple of) the driving frequency in the limit of vanishing energy (or as X — > 0) . In other 
words, expanding r]g^{X, r) as a Fourier series, we must have 

?7,^(X,r)- VA„(X)e™" with V|n|^l,A.(^) ^0. (43) 

X — *-0 

nGZ 

This implies that we have to investigate the first instability zone of the Hill's equation (|41|l . 
At a strictly zero energy {X = 0), this imposes A — for in this limit, Mlfl becomes 

Vq. + [1 - A] r;^^ = , 

(we again use the dot as a derivative with respect to r). The previous equation has a frequency equal to 1 for A = 
only. Thus, in a diagram (A, e), the Arnold tongue we are concerned with starts from the point (A = 0, e = 0). 

To evaluate the boundaries of this instability zone, we need to solve H41f) by expanding rjq^ (X, r) and A as power 
series in X: 

oo oo 

Tiq^ (X, r)=Y, viT^ {r)X^ ■ A = ^ A(™) X" . (44) 

m— m— 

As for the Mathieu equation, removing the secular terms from (|41|l gives rise to two critical curves lS.s{X) and /S.a{X) 
associated respectively with the periodic time reversal symmetric solution rjq^^^iX, r) = r/q^^s{X, — r) and the periodic 
time reversal antisymmetric solution r/g^.a(A', t) — r/g^^a(X, — r). 



3. First transition curve. 



Differentiating H4()|l with respect to t and dividing by X shows that 

^ (^^^ +f\X)V"{Qo{X,r)) (^f ) -0. (45) 

As Qq{X,t) is time reversal symmetric and equivalent to cos(t)X + 0{X) when X tends to zero (see (I15|) and 
(|22|l '). Qq{X,t)/ X is time reversal antisymmetric and satifies (|43|l . We then deduce that, for all X, Qq{X,t)/ X is 
an eigensolution of H41|) for the eigenvalue A = 0. In other words, riq^,a{X,T) = Qq{X,t)/ X is the eigenfunction 
corresponding to the eigenvalue A = whatever the energy (or X) is. We have then shown that the boundary Aa{X) 
of the instability zone we are concerned with is the e-axis of the A — e plane. 



4. Second transition curve for a partial isochronism of order n. 



To calculate the second curve, As{X) related to the time reversal symmetric solution 77,^^5 (X,t) of (|1T)) . we differ- 
entiate (|4U|) with respect to X: 



^ [ dQo 
9t2 1^ dX 



+ f\X)V"{Qo) 



fdQo\ , df^X)-, 



V dX 



dX 



V (Qo) = 0. 



(46) 



This equation is not of the Hill type because of its last term which is not linear in dQo/dX. 

Suppose now that the potential V is isochronous up to order n. By definition, we have T'^{X) = 1 + 2T2„ AT^" + 
o(A2"). That is ^ = 4nf2«A2»-i + o(A2»-1) and V' (Qo) = XQ^^^t) + o(X2). The last term is then of order 
A2«. 

This means that, up to order A^""^, equation H46|l reads 



\- 


(dQo\ 




[dXj 



dQo 
dX 



= 0, 0<fc<2n-l 



(47) 
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(Here and hereafter, the symbol [f{X)]k will denote the terms of order X'^ of f{X), [f{X)]k = hZx^(^))- Si'^ce 
dQo/dX = cos(r) + 0{X) as X ^ 0, it satisfies and is moreover time reversal symmetric. Now, because of (071, 
up to order X^"^^, dQo/dX represents the time reversal symmetric eigenfunction r]q^^g{X, r) of (|1T)| for an eigenvalue 
As which, up to this order, is equal to 0. This means that 



A,(X)= ^ A('")X™ and vi':^^) ^ [v,UX,r)l = 

rn—2n 

Let lis write now 



dX 



Q<k<2n-l. 



5Qo 
dX 



-C(r), 



(48) 



(49) 



2n 



where C,{t) is a function to be determined. Taking into account H48() and (I47II . we are able to calculate the terms of 
order X^" of ||1T|) and obtain 



^^^1^ + f\X) {V" iQoiX, r)) - a} r;,, {X, r) 







2n 



^ C(t) + C(r) - A(2") QW(r) - 4nf2„ Q^^^r) - 



(50) 



As Q^q\t) — cos(t), removing the secular terms from this last equation requires A 
of the second branch of the instability zone 

As{X) = -4nf'2„X2" + 0(X2") . 



(2n) _ 



-AnT2n- Hence, the equation 



(51) 



5. Bifurcation energy. 



Coming back to the definitions, A ~ (ljq q— o;^^ g)/i;2, X = (2e/u2)^/^, f'^{X) — V2/ujI andcjg = o^q g+7o,„ £"+o(e") 
(wg Q — U2), we can use (|51|) to obtain the leading order expression for the energy £g'^ at which a nonlinear in-phase 
mode isochronous up to order n bifurcates: 

(52) 




Let us first notice that the existence of a critical energy at which the in-phase mode undergoes a bifurcation implies 
that 7o,n is negative which, in turn, requires the corresponding frequency to decrease with increasing energy. This 
condition, already known for non isochronous potentials (or, in our terminology, isochronous up to order 1) [l3|, still 
holds in case of partial isochronism. 

As already mentioned in the introduction, expression (|52|l shows a very deep relation between the low-energy 
behaviour of the frequency of the in-phasc mode and its bifurcation energy. We will see in the next section that the 
same relation holds in fact for the out-of-phase mode. As n may now take on arbitrary positive integer values, H52|) 
treats the case of all possible analytic potentials with nonzero harmonic frequency {v2 > 0) . Its validity is however not 
restricted to such potentials and we remark that it still holds for the class of non analytic potentials studied recently 
by M. Kastner in relation with the possible existence of energy thresholds for discrete breathers [23ll2^. From this 
point of view, this relation seems quite general. 



6. Examples. 

The coefficients 70,0 — —2v2{2 / V2)"T2n are easily obtained from the expression (|A6|I for the period T. Indeed, eq. 
(|A6I) shows that T2„ = cr2n+i (2n -I- l)!!/(2n)!!. Then, zeroing the a2k+i, fc G {1, . . . , n — 1} as explained in Appendix 
IXIand re- introducing the corresponding constraints in (T2n-i-i, allows us to derive the expression for T2„ for a potential 
isochronous up to order n. 
For n = I (non isochronous potential), 

- 3 3 / 5 2 1 A 3vi 5v1 
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n=1 
n=2 
n=3 




FIG. 1: Schematic representation of the bifurcation energy 14111 as a function of the inverse squared number of oscillators A'^ 
for various degrees of isochronism of the in-phase mode (n =1,2 and 3). For a given n, the instability zone lies between the 
energy-axis (first transition curve, = 0) and the curve labelled by n = i, i =1,2 or 3. Below (above) the latter, the in-phase 

mode is stable (unstable). Right on a transition curve the solution is periodic. We see on this graph how the instability zone 
shrinks as n increases. In the purely isochronous case, the second transition curve merges with the first one on the energy-axis 
and the corresponding instability zone disappears. 



then, 



which is exactly the expression (3.13) of [T^. 

For n = 2, 



(1) _ 1202^1 sin^(^) 



~ 10 , ~ 15 f 1 28 4 ? \ 15 ^7 28 vj 1 ve 



Hence, 



6 81 ^ 15 V ' ■ V2 Vl5 ^2 81 6 1)2 



_(2) _ / 54z;|02 ,. /TT 



378w3^^5w| - 280w;f - 135t;6w| / V TV 



V2 sin - . (54) 



7. Pure isochronism 

Let us suppose now that the in-phase mode is purely isochronous. Then, T{X) — 1. Differentiating (|40|l with 
respect to X as we have done in the previous section gives rise to 

^'^"^ V"mi'-^^=0. (55) 



\dX V dX 
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At variance with eq. H46(l , this equation is of the Hill type and represents the time reversal symmetric eigensolution 
of the instability zone under investigation. Its eigenvalue is = for all X, i.e. whatever the energy is. We have thus 
found that, in case of pure isochronism of the onsite potential, the boundaries Aq(X) and As(X) of the instability 
zone merge and make it disappear. The merging of two transition curves is a phenomenon known in stability theory 
as a coexistence. It occurs for instance in Ince's or Lame's equations and general conditions for its appcarence are 
given in [2^1 12^. We notice incidentally that for an "isotonic" potential, V{x) — uP'lx + 1 — 1/ {x + l))^/8, which 
is isochronous (see for example the variational equation H41(l reduces to an Ince's equation. This allows for an 

explicit verification of the general result stated above. 

As a consequence of the disappearance of the instability region, the in-phase mode never undergoes a tangent 
bifurcation. Remark that, at variance with the results obtained for a partial isochronism, this result is non perturbative. 
Indeed, eq. and ifS^ are exact and their respective eigensolutions rjq^^aiX, r) = Qo/X and r]q^^s{X, r) = dQo/dX 
as well. 

To conclude, let us notice that expression (|52() shows how the instability zone shrinks as n tends to infinity (see 
also Fig. nj. Nevertheless, as it is valid in the limit of vanishing energies only it couldn't have been used to prove the 
result above. 



C. Out-of-phase mode 



Stability analysis. 



Evaluated along orbit II, the Jacobian matrix f^7\i now reads 



DF^kiQii) 



1 r 

2 

1 r 



k 



V'iQo + Q.) + r'(Qo - 0.) + 4sin^ - {IF"(2Q,) + W\-2Q^)} 



r'(Qo + Q.) - ^"(Qo - Q.) - 2zsinfc{W^"(2Q^) - W\-2Q^)} 



Sk,q 
Sk,q+Tr : 



(56) 



and the corresponding dynamics. 



1 r 
2 

1 r 

2 



V^" (Qo + Q.) + 1^" (Qo - Q.) + 4 sin^ I {I^" (2Q,) + W^" (-2Q,)} 
1^" (Qo + Q.)~ V" {Qo - Q.) + 2i sin q {ly" (2Q,) - W" (-2Q,)} 



Vq 

Vq+TT = . 



(57) 



To evaluate the possible bifurcating energies in this case, we will proceed along the same lines as for the in-phase 
mode. Notice, however, that eqs. H57|l do not decouple in the present case and that, due to the very last term, they 
are complex (real and imaginary parts are not decoupled) . To deal with this first difficulty, we will write the system 
(|57|l in terms of real and imaginary parts. 

Similar to the in-phase mode, the out-of-phase mode will eventually undergo a first tangent bifurcation via the 
perturbation rjq^ whose frequency is the closest to oJtt.q, that is, for Qc ^ ~ 2tt/N. As N tends to infinity, q — 2tt/N 
plays the role of the small parameter in the variational equations. But at variance with the in-phase variational 
equation where A oc sin^(7r/A'^) was the unique small parameter, eqs. (|57|l possess two small parameters through 
sin^(((7c + 7i')/2) ~ [tt/NY and sin((7c) ^ ir/N . Notice that these two parameters are not of the same order. 

In what follows, for the sake of clarity, we simplify further the notations by using rjq^ = r^-|-ij^ and rjq^+Tr = fo + ijo 
(with i = V— 1) rather than the lengthy values for qc and qc + vr. We work with dimensionless equations and with 
the same notations as in section Hi CI We take S = 4sin^(7r/A^) as small parameter. The quantities Tq, jo, ''tt, j-n-i Qo, 
Qtt and T depend on the amplitude y defined in III CI As the previous quantities, the small parameter S has to be 
expanded as a power series in y. Moreover, we define 



^2(r) 



V'' {Qo + Q.) + V^" (Qo - Q.) + (4 - <5) {W' (2g,) + ly" (-2Q,)} 
V''{Qo + Q-.)-V''{Qo-Q^) 



2 
2 

fv^(i- V4) [vF"(2g.) - w^"(-2g,) 

^v" (Qo + Q.) + f " (go - g.) + 4(5 {ly" (2g,) + (-2g,)} 



(58) 
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System (|57|) now reads (from now on, dots denote r derivative), 

' + Ti{6, T)r^ + T2(T)ro - %{5, T)jo = 0, 

j^+Ti{6,T)j^+T2{T)jo+T3{S,T)ro = 0, 

(59) 

f'o + 74(5, r)ro + T2(r)r^ T)j.^ = 0, 

^ jo + Ti(5,T)jo +T2{T)i.„ ~%{5,T)rT, = 0. 

As the energy (amplitude y) tends to zero, and Qo tend to zero as well and from (|58|l . H31() and H32(l . we obtain 
Ti{5^t) — > 1 — 5/3i/4, 7^((5, r) ai+ SPi/A and 7^(5, r) and T3((5, r) both tend to 0. As (5 is a small parameter, the 
only way to obtain a perturbation whose frequency is exactly 1 (as for the mode itself) is to require (5 = for it gives 
rise to 7i = 1 which is precisely the desired frequency for the perturbation (rjr, j^). Regarding {ro,jo), as 74 = ai < 1 
in this limit, we have to require (rg = 0,jo = 0) in order zero in y. The origin of the instability zone is thus located 
at i5 = 0, e = in the parameter space. 



2. First transition curve. 

Now, proceeding as in the previous section, let us derive the equations obeyed by Qo = Qo/y and Q^r ~ Qir/y- 

Q, + ri(o,T)Q, + r2(T)Qo = 0, 

(60) 

Qo + T4(0,t)Qo + T2(t)Q. = 0. 

As the out-of-phase mode has been choosen time reversal symmetric, Qo and Qt^ are time reversal antisymmetric. 
Moreover, due to the asymptotic relations given above, they represent a motion of frequency 1 as y — > 0. Having the 
required properties, they form the antisymmetric eigensolution of the Hill's system (|59|l for a value of 5 equal to zero 
for all y (or energy). Therefore, (r^,j,r) oc Qjr and (rp, jo) Qo- We thus deduce that, in a diagram (5,e), one of the 
boundaries of the Arnold tongue is the e-axis itself which means that, similar to the result obtained for the in-phase 
mode, no bifurcation can be expected from this branch. 



3. Second transition curve for an isochronism of order n. 



The second branch of the instability zone is the curve associated with the time reversal symmetric eigenfunction. 
Differentiating H30I) with respect to y, we obtain 







< 


\ dy J 








\ dy J 



ri(0,T) 

^4(0, r) 



V dy 

V 9y 



+ ^2(r) 
+ ^2(r) 



dy 

f dQ^ 
V dy 



dy dr-^ 

dlnf2 d^Qo 
dy dr'^ 



= 0, 



= 0. 



(61) 



Again, the functions dQ-^/dy and dQo/dy have the right time symmetry and Fourier properties. Nevertheless, the 
system above is not a Hill's system due to the last terms. For an out-of-phase mode isochronous up to order n, 
T{y) = 1 + r2„j/2" + o{y'^"), and we can evaluate the order of the last terms in (|61|l . Using (|31ll and H32|) . we find 



4nT2„ cosTi/2" + o(y 



dlnT2 d^Qo 



0{y 



2n+l\ 



(62) 



dy OT'' dy ot'^ 

Given this last result, we will first assume, and then verify a posteriori, that the order of 5 in y is 2n, hence 
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We then obtain immediately 



7^(r) = 2a2 cos Ty + o(y), 
r3((5,r) = v4(2^/32Cosry"+i +o(t/"+i), 
r4(5,r) = r4(0,r) + ^<5(2")y2« + o(2/2").. 



(64) 



We deduce from these relations that system (|59|l verifies 

f f,+ri(0,rK+r2(r)ro = 0(2/"), 

jV + ri(o,T)j, + r2(r)jo - o(y"), 

fo+T4(0,r)ro+T2(r)r^ = 0(2/"), 

I JO + r4 (0, T)io + (t)j„ = 0(2/") . 



(65) 



We remark now that the solution {riq^{T), 77g^+7r(T)) of (|57|1 is defined up to a global phase only (i.e. if {rjq^ (t), ?7qe+7r(T)) 
is a solution of lf57|l . (e^'^ry^^ (r), e*"^?7q^+7r(T)), (ji S M is also a solution). Using this freedom, we solve 0^ by requiring 











\^Qo^ 




dy 


k 


^0 — 


dy 



=0 ; = , l<k<n, 



(66) 



without loss of generality (the brackets have the same meaning as in the previous section). Reinstating in H59|) and 
using we find 



dy 



' '0 



dQo 
dy 



1< k <2n-l. 



(67) 



As in the last section, we now solve the first and the third equation of l|59f) at order 2n by introducing two functions 
Cti{t) and CoIt) according to 



' 7T 

^From H61|l and (|62|l we have 

9' fdQ^ 



dy 



\ dy 



71(0,' 



\ dy 



+ ^2(r) 



dy 
dQo 



+ Co = Co • 



(68) 



V dy 



and from ((231 and 



dy 

+ 'r2{r) 

/3i 



4?iT2„ cos r = 0, 



2n 



(69) 



= 0. 



2n 



Therefore, 



[f,+ri(0,r)r,+r2(r)ro]2„-^<5(2")cosr = 0, 
[fo+r4(0,r)ro+r2(r)r^]2„ = 0. 



C. + C^ - (4^5^'") + 4nf2„) COST = 0, 
Co +aiCo = 0. 



(70) 



(71) 



Removing the secular term from the first equation leads eventually to the equation of the second transition curve 

^ = -i^?/^" + o(2/^"). (72) 



/3i 



15 



Discarding the solution of the homogenous equation l^lj gives = 0. The second equation is satisfied by setting 



„(2n) 



and r. 



(2n) 



2n 



dQo 
dy 



Remark that these 



2n 



Co = as well (because its frequency ai ^ 1). Hence, ' 

results are correct provided we are able to justify a posteriori the initial assumption regarding the small parameter 
d. This amounts to proving that the remaining equations for and Jq in (|59|l never develop any secular terms up to 



order 2n. This is done in appendix |B| by showing that these equations verify the adequate solvability conditions. 



4- Bifurcation energy. 

Using the leading order relation between the amplitude and the energy density, — 2e/ajJ Q+o(e^), ui'^ q — U2+402, 
together with the relation 

iol = cjIo + 7.,ne" + o(£") , (73) 
and noting that 6 = 4sin^(7r/iV), f2ny^" = -7^.„e"/2ti;^ „ + o(£"), we finally obtain from {721 



/,,2 ,,2 \ . 2/^N \ 1/" 




(74) 



which is similar to the expression obtained for the in-phase mode. We notice this time that the existence of a 
critical energy at which the out-of-phase mode undergoes a bifurcation implies that 70. „ is positive. This requires 
the corresponding frequency to increase with the energy. This condition, already known for non isochronous motions 
p^ . still holds in case of partial isochronism. 

As seen in section III C 21 the period of the motion is more conveniently expressed as a series in the square of the 
amplitude, y^, rather than as a series in the energy e. For this reason, we also provide the bifurcation energy in terms 
of the coefficients T^n defined in lpTT]l , 



2 I nujl^gf^n j 



(75) 



5. Examples. 

We are now in a position to give an explicit expression of ei"'' for an out-of-phase mode isochronous up to order n 
in terms of the parameters Vi and 4>i defining V{x) and W{x) (see Q and (|28|l '). This amounts simply to finding the 
leading order of the low-energy (or low-amplitude) behaviour of the corresponding frequency as explained in section 
III C 21 (see also section IlII B|l . 

For a non isochronous out-of-phase mode (n = 1), from (|32|l and (|75|) . we find 



(1) _ A{v2+^2)(t>2 2 

3(^4 + 1604) + 3,, + r60. - it 



(^) (™) 



Let us notice first that, when f 3 = this expression reduces to formula (3.20) of (IJ obtained in the special case of a 
symmetric onsite potential V{x). The correction introduced by the asymmetry of V{x) (i.e. the term proportional to 
vl in the denominator of (|76|) ') has the interesting feature to be always negative. Therefore, the following inequality 



2 3112(^4 + 16(/'4)(3W2 + 1602) 



^ 2(5«2+3202) ^^^^ 

has to be satisfied for the out-of-phase mode to undergo a tangent bifurcation. This corresponds, as we have seen in 
the section above, to requiring that the frequency increases with the energy. 

Another interesting result easily drawn from (|76|l concerns the case of partially isochronous onsite potentials. It is 
found in this case that a certain amount of nonlinearity (^4) in the interaction potential is needed in order to ensure 
a bifurcation of the out-of-phase mode. Indeed, the relation between the first coefficients of the Taylor expansion of 
V{x) is lOvl = 9v2 V4 in this case. The denominator of H76|l is then positive provided 

'^4 > 7o rTFT' 

5 3V2 + 1602 
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So that, in a chain of harmonically coupled partially isochronous oscillators, no discrete breather (if any) stems from 
the tangent bifurcation of the out-of-phase mode. At the same time we can conclude, that breathers appear for fully 
isochronous harmonic oscillators (^4 = 0) when coupled anharmonically (04 > 0). Indeed, it has been recently proved 
that breathers exist and can be continued from zero anharmonic coupling for harmonic oscillators j29j . 
For a first degree of isochronism (n = 2), we obtain from H75I) 



2{V2 + 4(/)2)T4 



in(^), (79) 



where 



~ _ (15a^ + 96- 56ai)a| 9(32 + 5a? - 24ai)a3ai 3(5ai - 12)a4a2 5 

- 24a?(ai-4)2 IRR^^P ^ 8ai(ai-4) 16 + 

with an — Vn+i/{v2 + 4^2) and /3„ — 2"+^(/i„+i/(u2 + 4(/)2). Notice that, to enforce a partial isochronism of the 
out-of-phase mode up to order n = 2, we have used the first of the relations (|35|l {(3^ = f({aj})) to derive (jSUJ. 
Again, a tangent bifurcation of this mode will take place provided the inequality r4 < is fulfilled. 



6. Pure isochronism 



Similar to the pure isochronism of the in-phase mode, the pure isochronism of the out-of-phase mode leads to a 
merging of the two transition curves which makes the instability zone disappear. This is clear from the fact that, as 
f{y) = 1, dQ^/dy and dQo/dy obey 













(dQo 


^ 9r2 


[ dy 



ri(0,r) 



dy 

dQo 
dy 



^2(r) 



72 (r) 



dy 
dy 



0, 







(81) 



and are thus exact solutions of the variational equations H59|) for g = (i.e. (r,r) Jw) oc dQ-j^jdy and (rg, jo) BQ^/dy). 
As a consequence, if the out-of-phase mode is isochronous it doesn't undergo a tangent bifurcation. 



V. SUMMARY AND DISCUSSION 



We have examined so far the close link existing between the possible tangent bifurcation of band edge modes and 
the low-energy behaviour of their frequency. We have introduced the concept of partial isochronism of order n for 
these modes through the relation tj^ = ^2 ^ _(_ 7q_„e" -|- o(e") (g = or tt). By performing a linear stability analysis in 
the limit of small oscillations, we have derived a simple and general expression for the leading order of their bifurcation 
energy in terms of 7g,n. We have shown that the calculation of 70, « (in-phase mode) simply requires to invert a series 
and thus reduces to a pure algebraic problem. The inversion of this series is easily implemented with the help of any 
software able to perform formal calculations. The coefficient ^Ti.n related to the out-of-phase mode may be derived 
by means of a Lindstedt-Poincare expansion for the motion. 

In addition, we have proved by means of an exact linear stability analysis that, fully isochronous band edge modes 
(i.e. with an energy-independent frequency) never undergo a tangent bifurcation. At variance with the results quoted 
above, this one is non perturbative. 

In order to discuss the implications of these results on discrete breathers, we shall now assume that, at least some 
of them stem from the tangent bifurcation of the band edge modes investigated so far. To our knowledge, for a Klein- 
Gordon lattice with smooth (say C°°) linearizable onsite and interaction potentials (a;o,o and lUj^.o both nonzero) no 
general result exists which proves such an assumption. However, with the additional assumption that, in an infinite 
lattice, the breather amplitude (measured at the level of its largest oscillation) can be lowered to arbitrary small 
values |l3l Il5l|. we come to the conclusion that, in this limit, the breather frequency tends to an edge of the phonon 
band and merges with the corresponding mode. Indeed, for small amplitudes, the motion enters a quasi-linear regime 
and then approaches some phonon mode. But in the same time, the breather family (parametrised either by its 
amplitude, its frequency or its energy) has to lie outside the phonon band to avoid a resonance which would lead to 
its disappearance. And this is possible by approaching an edge of the phonon band only. Notice however that one 
can find systems with breather families which do not possess any small amplitude limit (see e.g. H^). Such breather 
families are then not related with the instabilities discussed above. 
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A. Energy thresholds for discrete breathers 



We now turn to an important implication of our results on the possible existence of energy thresholds for families of 
discrete breathers bifurcating tangently from band edge modes. It has alreadybeen noticed in many places that such a 
nonzero activation energy for discrete breathers is of practical relevance [Tsl 12^ as it surely affects their experimental 
detection and presumably their contribution to thermodynamical properties of lattices (see for example |28l | for some 
work in this direction). 

Let us consider a Hamiltonian lattice of N coupled oscillators described by with families of periodic orbits 
parametrised by their amplitude. Let us follow the periodic orbit corresponding to a band edge mode as its amplitude 
increases from zero. As we have seen in section Hvl under certain conditions, at a finite critical amplitude (or energy 
density) this mode will become unstable and bifurcate tangently. A family of discrete breathers emerging from this 
bifurcation is, right at this point, identical to the mode from which it stems and therefore has the same energy. 
Increasing the amplitude further and following this new orbit, we obtain its energy as a function of its amplitude. 
It is worth noticing immediately that in finite systems N < oo, discrete breathers arising in this way exist above a 
certain energy threshold only. Indeed, due to their finite amplitude their energy is surely nonzero. The question is as 
whether this threshold persists in the thermodynamic limit, N —^ oo. 

We then turn to the determination of the (total) bifurcation energy E^q'^ of a BEM partially isochronous up to 
order n, in the thermodynamic limit. According to (|52|l and l|74|l . the bifurcation energy density for the two possible 
modes q = or q = n can be cast into the general form 



2 2 



1/n 



where the coefficient 7g_„ is defined by the low-energy behaviour of the BEM frequency 

i^l = ujI, + 7,,„£" + o(£") , (83) 

and where the wave number q^ is the closest to q { q^ = 2t:/N if g = and q^ ^ n — 2t:/N if q = n). In any case, 
N-^. We then find that 



/ ,2 , ,2 

^9,0 ~ ^q^,0 



=iVe(«) ^iVi-t (N^oo). (84) 

This means that if the BEM is not isochronous, (n = 1), its total bifurcation energy vanishes as the lattice becomes 
infinite and so does the breather energy in this limit. No energy threshold exists in this case as already mentioned in 
|l3|. However, as soon as the band edge mode bears some degree of isochronism, (n > 1), its total bifurcation energy 
either converges to a finite value (n = 2) or simply diverges {n > 2) and energy thresholds are thus expected. 

We note incidentally that expression (|84ll . although valid for a one-dimensional chain, bears some striking resem- 
blance with its multi-dimensional counterpart in the non isochronous case which reads Eq'^^ ^ N^~i where d is the 
dimension of the lattice [l^ [l5l | . Combination of both isochronism (n) and dimensionality [d] leads immediately to 
the conlusion that the total bifurcation energy of a BEM scales like 

(TV^oo). (85) 

Energy thresholds for discrete breather families bifurcating tangently from BEMs are thus expected as soon as one of 
the positive integers n or d is strictly greater than one. 

As already mentioned in llVB 5l the general expressi on 118211 still holds in case the analyticity of the onsite potential 
is relaxed. Recent results obtained by Kastner in [23l l24| prove for example that, for (ID) onsite potentials of the 
form V{x) — /2 + \xY /r + o(|a;|''), where r is any real number greater than 2, H82|l and 1)83(1 are still valid and 
that n = r/2 — 1. This way, the exponent n governing the low-energy behaviour of the frequency (|83|l can be tuned 
continuously over the whole range of positive real numbers and may counterbalance the effect due to the dimensionality 
d. This nonintegral partial isochronism opens up the possibility to satisfy the inequality obtained from H85|l 



nd<2, n e K+ , d e N+ , (86) 

ensuring the absence of energy threshold in dimension d, even in two- or three-dimensional systems. 

Finally, we also mention the existence of energy thresholds for breathers in one-dimensional systems with alge- 
braically decaying long range interactions |3Clj | . We can then identify three lattice properties which lead to appearance 
of breather energy thresholds - dimensionality, interaction range, and (partial) isochronism. 
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B. Energy thresholds for discrete breathers revisited 

1. DNLS equations for the slow modulations of BEMs 



To conclude this paper, we propose to revisit certain results of the previous section by deriving a discrete nonlinear 
Schrodinger equation (DNLS) for the slow modulations of small-amplitude partially isochronous BEMs. We shall do 
it with the help of a method based on their nonlinear dispersion relation (see for example |3ll is^ . ) which renders 
its derivation almost straightforward. We would like to mention that this result, obtained by a somehow heuristic 
method, is also confirmed by a more rigorous multiple-scale analysis at least for the first orders of isochronism. 

The method we will use to derive the DNLS equation is not properly speaking a multiple-scale expansion for 
the motions xi{t) which would read xi(t) = X]j>i M"'^j./(^Oj ^ij ■ ■ ■ ) where /x <C 1 is a small parameter and where the 
different timescales are given by i,„ — (see for example 0) US HESS)- But it represents a similar approximation 
and has the advantage to be more transparent. It is based on the nonlinear dispersion relation obeyed by a BEM 
isochronous up to order n. 

We first note that in the harmonic limit, plane waves Xq-i{t) cx e^^'^^~^''-°^\ where lj^ q is given in (|10|) . are solutions 
of the linearised equations of motion J^. Looking for slow modulations of the q ~ and q ~ n modes, we write them 



Xg.i (t) = Aie^^^'-""'"*) + c.c. -I- oiii) (87) 

where "c.c." stands for "complex conjugate", q E {0,tt} and the small amplitude Ai = fiipi. Here is some small 
parameter. 

The function tpi is slowly varying in space and time. It is found typically from a multiple-scale analysis that, for 
non isochronous modes, it assumes the form ipi = tpidi^t, ^l) (see for example For modes isochronous up to 

order n, a similar analysis yields ipi = ipi{^'^"t, ^J,"l). This is easily understood from the nonlinear dispersion relation 
Lu'^ ~ Q -|- 7g,n£"- For a slowly modulated plane wave to exist, it is known that nonlinearity has to compensate 
the dispersion of the wave packet. Now, the nonlinear term in the previous relation is, in leading order, proportional 
to e" ^ /i^". Therefore, the nonlinear correction to the frequency introduces the natural long-time scale On 
another hand, the dispersion of the wave packet arises from a (spatial) Laplacian. If fx^l is the space scale, dispersion 
introduces a correction of order /i^^. This correction is able to compensate the nonlinear effects only if p — n, hence 
the scaling of tpi. 

Now, in leading order (harmonic regime), the energy density of a BEM reads = 2w^o|Ap, q e {0,7r} (the BEM 
amplitude A is constant so that it doesn't depend on /). Re- instating in the nonlinear dispersion and expanding the 
latter around q = or g = tt yields 

t^o-^o.o^ — [l-cos(g)] + -!^(2^2^)"|A|2" and u;^-^^^, ^ [l-cosiq-7r)] + ^{2u;lJ^\A\'^ . (88) 

As explained in 31, wc can now use the dispersion relations above to find the equations governing the envelope 

of the slowly varying amplitude A. This merely amounts to replacing [ujk — (^k,o) by idt, (q — k) by idi (with fc = 0, tt), 
A by the slowly varying amplitude Ai and finally to let it act on the amplitude Ai. We then obtain, 

idtAi = --^{Ai+i+Ai-i-2Ai) + ^{2ul^r\Ai\''''Ai, (in -phase) (89) 
idtAi = J^{Ai+^+Ai^^-2Ai) + ^{2ujlXWAi, (out - of - phase) . (90) 

Several remarks are in order at this stage. First, the above derivation is not rigorous. Nevertheless, it makes the 
role of isochronism through the nonlinear dispersion relation transparent and allows us to derive the DNLS equations 
for the modulation of the BEMs very easily. Second, we remark that for non isochronous potentials [n = 1), H89|) 
is exactly the equation (6) derived by Kivshar in |33| (see also 39J, eq. (5)). In this respect, (|89|l and H9U|) are 
generalisations to higher order of isochronism. We would like to mention as well that, for n = 1, (|89() and (|90|l have 
been obtained by a rigorous multiple-scale expansion in '37] (up to the fact that the discrete spatial derivative is 
replaced by a continuous one). 

Equations (|89|) and H9UI) provide now a physical motivation for studying DNLS with higher order nonlinearities. 
Indeed, the nonlinear exponent is shown to be directly related to the isochronism of the BEM under consideration. 
We notice finally that, (|89|l and H9U|) describe both the BEMs and their modulations and that, for this reason they 
make it possible to study the bifurcation process which leads from the plane wave to a breather (within the degree of 
approximation inherent to their derivation). 
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2. Rederivation of discrete breather energy thresholds 



Now, the condition for the stability of a plane wave for the cubic DNLS equation on a periodic lattice has been 
derived by Carr and Eilbeck in Studying an equation of the type iAi + x(^/+i + M-i — 2^;) + A|A;pA/, with 
^i+AT = Ai, N being the number of lattice sites and N = l^/P = these authors show that the plane wave 
Ai = e^'^^A is stable if and only if A < 2iVx sin^(7r/Af). The generalisation of this result to arbitrary nonlinearity 
(A|A/P"A/) and norm J\f is straightforward and reads 

l/n 

(91) 



TV - 



2X 2, 



, sin^(7r/iV) 



For a plane wave Af / N — \A\'^ and according to the mode q under consideration, the amplitude and the energy density 
are related via Eg — 2a;^Q|y4qp. Moreover, in H89|) and (|90|l . the ratio of the coupling Xq and the nonlinear parameter 
Xq is Xq/^q = ±02/(2Wg o)"7q.n , "+" for q = TT and "— " for q — 0. Then, eventually, we find the energies at which 
the in-phase and out-of-phase modes become respectively unstable 



which are precisely the expressions obtained in the previous sections. Notice that, these energies exist provided 
70,71 < and j^^n > respectively. At the level of H89|l and (|90|l . these conditions imply Xq\ > Oj (9 = 0, tt), which 
is precisely the condition for these equations to support bright type breathers. 

Our concluding remark is that, for nonlinearities higher than the usual cubic one, the DNLS equation is known to 
possess energy thresholds for discrete breathers (see for example [TEl IT^ llllli^ ). This allows us to corroborate our 
previous prediction that discrete breathers steming from the tangent bifurcation of partially isochronous band edge 
modes appear above certain energy thresholds only. Equations (|S^ and (pn|l should provide a way to estimate them. 



VI. CONCLUSIONS 



In this paper we have shown that, in network of coupled oscillators described by Hamiltonian (Q, the way band edge 
modes (BEMs) possibly undergo a tangent bifurcation crucially relies on the low-energy behaviour of their frequency. 
This behaviour is obtained by expanding the frequency as a power series in the energy (see (jHSI) as explained in 
section III B 21 and in appendix 1X1 for the in-phase mode or in section III C 21 for the out-of-phase mode. 

The energy at which such modes may bifurcate tangently is given in (|52|1 for an in-phase mode and in (|74|) for an 
out-of-phase mode. In these expressions, n is the degree of isochronism of the mode under consideration and 7g,,i, the 
first nonzero nonlinear coefficient of the low-energy expansion of the frequency as defined in (|83|l . 

As proved in sections IIV B 71 and IIV C 61 purely isochronous modes (i.e. with an energy-independent frequency) 
never undergo a tangent bifurcation. This result is non perturbative. Now, in a network of coupled isochronous 
onsite potentials, the in-phase mode is isochronous and never bifurcates tangently according to the previous theorem. 
But the out-of-phase mode, which is generically non isochronous, may bifurcate as soon as the nonlinearity of the 
interaction potential is strong enough as expressed in H78() . 

Bifurcation energies obtained in (|52|l and (|74|l are used in section IV AI to derive a condition for the occurrence of 
energy thresholds for discrete breathers bifurcating tangently from BEMs. For modes isochronous up to order n and 
in dimension d, this reads nd > 2. For analytic potentials {n positive integer) energy thresholds exist whenever n > 2 
(even in ID) or when d > 2. For non-analytic potentials, n becomes a positive real number which opens up the 
possibility of absence of energy threshold even in 2- or 3D. 

In section IV B II we derived two DNLS equations respectively related to the slow modulations of small-amplitude 
partially isochronous in-phase (eq. (|HH) and out-of-phase (eq. (|90|) ') modes. These equations are the generalisation 
to higher degree of isochronism of the usual cubic DNLS equation. Their nonlinearity is proportional to the degree 
of isochronism of the corresponding mode As shown in section IVB 21 they allow for an easy rederivation of the 
bifurcation energies H52I) and H74|l . 
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APPENDIX A: ENERGY EXPANSION FOR THE PERIOD IN A ID POTENTIAL 



We derive in this appendix a simple method to compute the period of osciUation around the minimum of a potential 
V{x) as a power series in the energy E. This method is, to our opinion, simpler than other methods which have been 
developed earlier 0, ^|. It yields the coefficients of this power series in terms of the coefficients of the Taylor 
expansion of V{x) around its minimum. This is used to construct a potential partially isochronous up to order n, 
which merely amounts to canceling out all the coefficients of the power series in E up to order n — 1. 

We assume that V{x) has the following properties: V{0) = 0,V'{0) = 0, y(0) = ^ 0. Moreover, we require 
that V'{x)/x > on an interval / centered around the minimum. In this interval, V{x) can be expanded as a power 
series in x around 0, 

1 °° 1 

Vix) = -uj^x^ +uj^Y -an-ix''. (Al) 

2 ^ n 

At a given energy E, the period T is given by 

r- r+'^^^ dx 

T{E) . x/2 / /\ (A2) 

Jx-(E) ^/E-V{x) 

where x±{E) are the two solutions of V{x) — E, {x+{E) > and x^{E) < 0). Let us make a change of variable 
defined by V{x) — ^uj'^X'^, with dX{x)/dx > for x E I. The last condition ensures that AT is a monotonic increasing 
function of x. Then (|A2|) yields, 



TiE)^V2 " -^^S±^dX = ^ A'^u]^^ (A3) 




where A(A) = x{X) — x{—X) with dx{X)/dX > 0. A' represents the derivative with respect to X. In order to 
compute the last expression of T{E), we need to invert the change of variable. Taking into account the condition 
dx{X)/dX > 0, we obtain 

-J oo 

V{x) = 2^^^^^ = ^ + H (^"^^ 

where the coefficients (t„ are functions of the am- Hence 

A'(A) = 2 + f2{2n + 1)(72„+iA2" j . (A5) 
Reinstating this expression into HA3|I and performing a term by term integration over u, we finally get 



TiE) = ^ 

UJ 



' (2n+l)!! /2^\" 

1 + a'2„+i 



(A6) 



(2n)!! \oj^ J 

The inversion HA4(I is easily done with the help of any software able to perform formal calculations. Using Maple [20| . 
we have computed the few first odd coefficients (T2n+i- 

5 2 1 

as = — an as 

18 2 4 ^ 

77 4 7 2 7 7 2 1 

^ 216"^ ~ 8"^"^ + 15"^"^ + 32"^ ~ 6"^ 

3 33 143 o 9 2 715 , 143 2 2 H 2 3 2431 g 33 , 1 

aj = -a,a, - -am + —a,a, + -a, - —a,as + —a,a, - -a,a, + -asa, + ^^a, - —a, - -a, 

To construct a potential partially isochronous up to order n we need to zero the coefficients (T2fc+i up to cr2n-i- This 
constraints the coefficients at to verify certain relations given below: 

Order n = 2 ; as = -TrOi\ 

9 

56 4 14 

Order n = 3 ; a^ = — — 02 + —0^40:2 



24 592 3 36 2 848 . 
Order n = 4 ; 017 = —aQa2 —0402 + ^^"4 + -^^2 
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APPENDIX B: SOLVABILITY CONDITION FOR SYSTEM JS^ 

We prove hereafter that, up to order 2?! in the amphtude y, system (|59|) never develops any secular term. 

• Up to order y^", we have shown in section HV C 31 that the real parts r,r and ro of the two modes ?/,r-27r/Ar and 
'>l2ii/N obey the system 

' f^+Ti(0,r)r^+T2(T)ro = 0, 
fo+T4(0,r)ro+r2(T)r^ = 0, 

The functions 7^(t) are 27r-periodic. 

• Up to the same order, the two periodic solutions of this system are given by 



(Bl) 



dy 

dQo 
V dy J 



and 



( 1^ \ 

y dr 

I dQo 
\ y dr / 



(B2) 



The first solution is time-symmetric ((^•') whereas the second is time-antisymmetric (^"^). 

• From order to n we can take the imaginary parts and jo of the modes equal to 0. From (n + 1) to 2n they 
obey, 



(B3) 



jV+Ti(0,r)j,+T2(T)jo+T3(5,T)r|,^^ = 0, 

'jo + %{0,T)jo+T2{T)j^~%i6,T)ri:^ = 0, 
where Ta is time-periodic and time-symmetric (it is a function of Q^^ which is time-symmetric itself). 

• The problem is to prove that frornj/""'"^^ to y^", (|B3|I never develops secular terms. We will use the solvability 
condition below (see for example [43). 

Theorem 1 (Fredholm alternative) Assume the subspace of T -periodic solutions of the homogeneous system 



-\x) ^ Ait)\x), 

A{t + T) — A(t), A G Mat(ri x n) is of dimension k > 1. Denote by |?7*-'''), I G {1, • • • , k}, k linearly independent 
T -periodic solutions of the adjoint system 



dt 



\v)=~A\t)\r,). 



For \f{t + T)) = \f(t)), the inhomogeneous system 

d 



dt 



\y)^A{t)\y) + \f{t)), 



has T -periodic solutions if and only if 



{r^^'^\fit))dt = [l e {!,■■■ ,k}). 



(B4) 







• In our case, we rewrite (IB3|I as a first order inhomogeneous system with \y) — col(r^,rQ,pT^,pQ) (column vector) 
and Ptt = Tt^,Pq = i-Q. We have 



A 



/ 10 

1, , 

-Ti -r^ I ^'^^ = 

V -Ta -% 





V %ri'^ I 



(B5) 
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The corresponding homogeneous adjoint system is 




Ti 
T2 % 
-10 
0-100 





(m\ 


) 


m 




\m 1 




\ ??4 / 



(B6) 



Its solutions are 



• Calculating the two solvability conditions (|B4p . we get 




(s),(a) 




27r 



27r 



(s),(a) 



and 



(B7) 



(B8) 



(^^^' I/) dt = / %{-r'^:\i'^ + r^'^Vi'')) dr = 



(B9) 



Indeed, to compute this last expression, we remember that is time-symmetric (series of cosine terms). The 
term in parenthesis being time-antisymmetric, the integrand is time-antisymmetric. Because it is 27r-periodic, 
the integral is zero. 

• Conclusion: The solvability conditions are fulfilled and consequently (|B3|I does not develop any secular terms 
up to order 2n. The way we have obtained the critical energy for the tangent bifurcation is thus correct. 
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